Spectroscopy system based on a single quantum cascade laser for simultaneous detection of CO and CO2
Wei Min1, 2, Ye Qing-Hao3, Kan Rui-Feng1, Chen Bing1, Yang Chen-Guang1, Xu Zhen-Yu1, Chen Xiang1, 2, Ruan Jun1, Fan Xue-Li1, Wang Wei1, Hu Mai1, Liu Jian-Guo1, †,
Anhui Institute of Optics and Fine Mechanics, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei 230031, China
University of Science and Technology of China, Hefei 230022, China
Shanghai Institute of Satellite Engineering, Shanghai 200050, China

 

† Corresponding author. E-mail: jgliu@aiofm.ac.cn

Project supported by the National Key Scientific Instrument and Equipment Development Project of China (Grnat No. 2014YQ060537), the National Basic Research Program of China (Grant No. 2013CB632803), the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant No. XDA05040102), and the National Natural Science Foundation of China (Grant No. 41405134).

Abstract
Abstract

A quantum cascade laser (QCL) based system for simultaneous detection of CO and CO2 is developed. The QCL can scan over two neighboring CO (2055.40 cm−1) and CO2 (2055.16 cm−1) lines with a single current scan. The wavelength modulation spectroscopy (f = 20 kHz) is utilized to enhance the signal-to-noise ratio. A white cell with an effective optical path length of 74 m is used. The calibration of the sensor is performed and minimum detection limits of 1.3 ppb (1 × 10−9) for CO and 0.44 ppm (1 × 10−6) for CO2 are achieved.

1. Introduction

Carbon monoxide (CO) and carbon dioxide (CO2) affect the heat balance of the atmosphere and participate in the atmospheric chemistry.[1] CO is an important atmospheric pollutant in the troposphere.[2] CO also serves as an indicator of urban pollution, and its anthropogenic source is mainly the incomplete combustion of hydrocarbon fuels.[3] Moreover, CO can react with O2 to produce CO2, therefore it is an indirect greenhouse gas (GHG). CO2 is the most important anthropogenic GHG. It is the main product of combustion processes and human activities. Its annual emissions had grown by about 80% between 1970 and 2004, from 21 gigatonnes (Gt) to 38 Gt, and accounted for 77% of the total anthropogenic GHG emissions in 2004.[4] Reliable and sensitive detection of atmospheric CO and CO2 simultaneously is necessary,[5,6] as these measurements can be used to determine their sources, leading to an improved understanding of the carbon cycle and climate change.

Tunable diode absorption spectroscopy (TDLAS) technology is ideally suited for trace gas detection in environment monitoring, industry control, and combustion diagnosis due to its advantages of in situ, high sensitivity, non-intrusive, and high selectivity.[711] Simultaneous detection of atmospheric CO and CO2 based on TDLAS has been reported by several research groups.[1215] Gabrysch et al.[12] used a near-infrared laser (1.578 μm) for CO and CO2 simultaneous detection and the sensitivity was increased by a factor of 15 compared to direct absorption, the minimum detection limit for CO and CO2 was about 10 ppm (100 torr, 1 m) and 53 ppm (8 torr, 1 m). Engelbrecht[13] reported a compact fiber-optic diode laser spectrometer for the two species measurement around 1580 nm, achieving a detection limit of 5.1 ppm CO and 9.1 ppm CO2 with a 1 m absorption path length at 80 hpa. Chen[14] developed a system for the detection of CO and CO2 with two near-infrared lasers (1563.62 nm and 1580.99 nm) and a white cell. A remote laser-based sensor for vehicle CO and CO2 emission detection was designed with a single laser (1.57 μm).[15] Those researches were all based on the near-infrared spectrum, however, the fundamental absorption lines of most trace gases are in the mid-infrared region.[16] Quantum cascade laser (QCL) is now a perfect mid-infrared light source[17] and it is more suitable for high-sensitive trace gas detection in combination with TDLAS and the wavelength modulation spectroscopy (WMS) technique.[1820] Recent developments of QCL for simultaneous atmospheric trace gases detection have been reported, such as simultaneous detection of atmospheric CH4 and N2O emission fluxes by QCLS (4.5 μm, 7.5 μm, and 7.9 μm),[21] and simultaneous atmospheric N2O and CO detection by a 4.5 μm QCL.[22] A QCL (7.73 μm) based sensor for simultaneous detection of atmospheric N2O, CH4, and H2O was developed with a Herriot cell at 100 torr.[23] However, using a single QCL for simultaneous measurement of atmospheric CO and CO2 has rarely been reported.

In this paper, we develop a mid-infrared laser sensor for simultaneous measurements of atmospheric CO and CO2 based on a single 4.85 μm continuous wave QCL (CW-QCL). The appropriate wavelength enables simultaneous multi-species measurements using a single laser, simplifying the sensor system and reducing its cost greatly. Two neighboring CO (2055.40 cm−1) and CO2 (2055.16 cm−1) absorption lines can be covered in a single QCL wavelength scan. A multiple pass absorption cell with an effective path length of 74 m and the WMS technique are employed to achieve high sensitivity.

2. Theory and absorption line selection
2.1. Basic theory of WMS

The absorption spectroscopy obeys the Beer–Lambert law. When adding a sinusoidal modulation of angular frequency ω to the laser, the time variations of the laser frequency ν and the laser intensity are expressed as

where is the laser center frequency, α is the frequency modulation amplitude depth, I0 is the average laser intensity, i1 and i2 are the linear and nonlinear amplitudes of the intensity modulation (normalized by I0), and ψ1 and ψ2 are the phase shifts between intensity and frequency modulation for linear and nonlinear terms, respectively.

The transmission coefficient τ can be described in terms of a Fourier series[24]

Combining Eq. (3) with the Beer–Lambert relation, the H terms can be defined as

where Sj(T) is the line strength of absorption transition j at temperature T, P is the total pressure, Xi is the mole fraction of absorption species i, L is the path length, and ϕj is the line-shape function of transition j. When the absorption is more than 0.1, the H terms can be written as

In general, a lock-in amplifier can be used to acquire the wavelength modulation 2 f signal by multiplying the detect signal with the reference sinusoid and applying low-pass filters to separate the 2 f component and eliminate other components outside the filter band. The 1 f normalized method is often used in WMS to eliminate the influence of light intensity variations.[25,26] In most cases, i2 is much smaller than 1, and we assume that the phase shift of the linear intensity modulation is π, then the 1 f-normalized 2 f amplitude can be simplified to

The absorption gas concentration X can be calculated according to the normalized 2 f signal I2 f−normalized if the gas temperature, pressure, and optical path are fixed. When the gas temperature changes, the measurement result can be corrected with the S(T)/S(T0) factor, where T and T0 represent the measurement and the reference temperatures, respectively.[27]

2.2. Absorption line selection

Two neighboring CO and CO2 absorption lines need to be selected for simultaneous detection. The absorption lines of CO and CO2 in the mid-infrared spectral region (2000–4000 cm−1) are shown in Fig. 1 based on the HITRAN database.[28] The absorptions of H2O and N2O in the same spectral region are also given, which are the potential interferences. For CO2 and CO, their strong absorption lines are within the spectral range of 2000–2400 cm−1, while N2O and H2O also have absorption lines within the range.

Fig. 1. The absorption line strengths of (a) H2O, (b) N2O, (c) CO, and (d) CO2 at 296 K in the mid-infrared spectra (HITRAN database).

A commercial CW-QCL laser (Alpes laser) covers a wavelength range of 2052–2058 cm−1 through tuning its temperature and current. The absorption spectra for 2% H2O, 380 ppm CO2, 200 ppb CO, 320 ppb N2O in air are simulated at 269 K and 250 torr to help to choose the suitable lines, as shown in Fig. 2. From Fig. 2, the two neighboring lines of CO at 2055.40 cm−1 (P21) and CO2 at 2055.16 cm−1 (P28(e)) can be selected, which are fully swept in a single wavelength scan at a fixed temperature (injection current between 140 mA and 165 mA). Their HITRAN parameters are shown in Table 1. Furthermore, the second harmonic signals of different mixture gases of CO and CO2 at 250 torr are simulated in Fig. 3 to verify that a relative lower pressure (250 torr) can be selected to avoid the spectral interference between CO2 and CO. From the simulated results, the amplitude of CO (200 ppb) varies ∼0.2% when CO2 changes from 380 ppm to 1000 ppm. Similarly, the amplitude of CO2 (380 ppm) varies ∼0.018% when CO changes from 200 ppb to 3 ppm.

Fig. 2. The simulated absorption spectra for H2O, N2O, CO2, and CO in air at 296 K and 250 torr based on HITRAN database.
Fig. 3. The simulated second harmonic signals for different mixtures of CO2 and CO.
Table 1.

Spectroscopic parameters of the selected CO2 transition at 2055.16 cm−1 and CO transition at 2055.40 cm−1 at 296 K.

.
3. Experimental system

The experimental single-QCL-based system is presented in Fig. 4. A continuous wave distributed feedback QCL (Alpes Laser, Switzerland) around 2055 cm−1 in a high heat load (HHL) package was used as the light source, with a beam-divergence angle less than 3 mrad. The temperature and the current were simultaneously controlled by a commercial laser controller (LDC 3724B, ILX Lightwave, USA). The mid-infrared light was combined with an aligned visible beam (λ = 650 nm) by a beam splitter (T > 90% at 4.85 μm, R > 85% at 650 nm) for aiding the optical alignment. The combined light beam went through two reflect mirrors, then passed into a white multi-pass cell. The cell consists of three spherical mirrors, with a base length of 1 m, which can provide an optical path length of 74 m. The leakage rate of the multi-pass white cell is 0.01 Pa·L/S. The output mid-infrared beam was focused onto an AC thermoelectric cooled mercury cadmium telluride detector (MCT) (PVI-2TE-10.6, Vigo Systems, Poland) using a parabolic mirror (f = 50.4 mm). The detected signal was acquired by a data acquisition (DAQ) card (M2i.4210 PCI 14 bit 20 Ms/s, SPECTRUM, Germany), and then input to a Labview-based lock-in amplifier to extract the 1 f-normalized 2 f signal.

Fig. 4. Schematic of the experimental apparatus.

The optical path outside the gas cell was about 50 cm, which was 148 times smaller than the inside optical path of 74 m, so the influence of CO and CO2 in ambient air can be neglected. It is better to make the outside part closed or purge with N2 particularly for CO2 detection, because it is easily influenced by human breath in the laboratory.

The two absorption lines were covered by a sweep range of about 0.6 cm−1 with a current scan signal of 200 Hz. The modulation depth was optimized, which affects the peak value of the 2 f signal. The modulation frequency was set at 20 kHz to reduce the low-frequency noise. Under this modulation frequency, the modulation depth varied with the modulation voltage directly, so an appropriate modulation voltage should be selected. The laboratory air flowed into the cell and the needle valves were adjusted to maintain the pressure of 250 torr. The peak value of the 1 f-normalized 2 f signal is plotted in Fig. 5. As shown in Fig. 5, both the 2 f peaks of CO and CO2 reached the maximum with a modulation voltage of about 120 mV, thus this modulation voltage was selected for CO and CO2 simultaneous detection in the following experiment. The detector signal, the 1 f-normalized 2 f signal, and the peak values of CO2 and CO are also given in Fig. 6. The asymmetry of the 2 f signal in Fig. 6 resulted from the intensity modulation. In Fig. 6, the 1 f-normalized 2 f signal of CO is different from that of CO2, this difference perhaps is due to the background effect. Also, the absorbance of CO in the laboratory air was much less than that of CO2, which would result in different background effects to their 2 f signals.

Fig. 5. The 2 f peaks of CO and CO2 in laboratory air at different modulation voltages.
Fig. 6. An example of detector signal and related 2 f signal of laboratory air.
4. Experimental results

Several standard gases (Guangming Research & Design Institute of Chemical Industry, with 2% relative uncertainty) with different mole fractions of CO (200 ppb, 500 ppb, 950 ppb, 1.7 ppm, 3.7 ppm) and CO2 (95 ppm, 257 ppm, 404 ppm, 702 ppm, 1001 ppm) were used for the sensor calibration, and the carrier gas was N2. During the calibration, these standard gases flowed into the gas cell with the concentration increased. The 2 f amplitude increased with the concentration of both gases. When each 2 f peak became stable, the signal was recorded for 20–30 min. The calibration processes of CO and CO2 are plotted in Figs. 7 and 8. After calibration, a continuous detection of relatively lower standard gas (950 ppb CO, 404 ppm CO2) was conducted to test the sensor reproducibility, as illustrated in Figs. 7(a) and 8(a).

Fig. 7. Calibration process (a) and calibration results (b) for CO.
Fig. 8. Calibration process (a) and calibration results (b) for CO2.

For CO, the linear fit coefficient of CO is greater than 0.998. For the second 950 ppb CO detection, the average 2 f amplitude in 30 min is 2.14, and the corresponding value is 2.12 for the first detection in the calibration, so the amplitude difference is 0.02. The related concentration difference between the two stages is 10 ppb (1%). The two differences are small, showing the good reproducibility. Similarly, for CO2, the linear fit coefficient is greater than 0.99. For two detections of the same concentration of the standard gases, the differences are 0.35 and 15 ppm (3.7%) for the amplitude and the concentration, respectively, which are both relatively small. All these results illustrate a good linear response and reproducibility of the sensor for the detection of CO and CO2.

In order to test the detection limit and the long-term precision of the experimental system, the laboratory ambient air was closed in the gas cell and continuously measured at 250 torr. An example of the signals is shown in Fig. 9. The continuous detection results of 3000 s are shown in Fig. 10 with 1-Hz sample rate. According to the measurement results, the averaged concentrations of CO and CO2 are about 1.07 ppm and 509.24 ppm, which are in the normal concentration ranges in atmosphere. During the detection period, the root mean square (RMS) noise is 8 ppb for CO (about 0.7% of the laboratory ambient air level) and 1.78 ppm for CO2 (about 0.35% of the laboratory ambient air level).

Fig. 9. An example of the related signals of ambient air: (a) detector absorption signal, (b) normalized 2 f and 1 f signals.
Fig. 10. Continuous measured results of CO (a) and CO2 (b) in laboratory ambient air for 3000 s.

In general, the Allan variance plot is used to seek the optimum integration time and the corresponding minimum detection limit for the system. The Allan variance of the simultaneous detection results is plotted in Fig. 11, which demonstrates that the optimum integration time is about 85 s for the simultaneous detection of CO and CO2. A minimum Allan deviation of 1.3 ppb for CO (∼ 1‰ of the ambient level) and 0.44 ppm for CO2 (∼0.9‰ of the ambient level) is achieved at the optimum integration time. The high detection limit shows that the system is appropriate for atmospheric monitoring of these two species. With a longer integration time, there is an increase in the Allan Deviation due to the instrumental drifts, including laser wavelength drift, optical drift, and electronic drift.

Fig. 11. Allan deviations of CO (a) and CO2 (b) as a function of integration time.
5. Conclusion

We have developed a single QCL based absorption spectroscopy system for simultaneous detection of atmospheric CO and CO2 with a 74 m multi-pass cell. Utilization of a single QCL to realize dual-species detection reduced the cost and simplified the experimental configuration. The two neighboring CO (2055.40 cm−1) and CO2 (2055.16 cm−1) absorption lines were covered within a single wavelength scan. The concentration calibration of the system was performed using standard gases with different concentrations, and the calibration results illustrated a rather good linear response for both CO (200 ppb–3.7 ppm) and CO2 (95 ppm–1001 ppm). The detection limit and long-term precision of this CO/CO2 sensor were examined by monitoring the laboratory room air for about 3000 s, and its Allan deviation plots showed that a detection limit of 1.3 ppb for CO and 0.44 ppm for CO2 can be achieved at 85 s integration time.

Reference
1Tang X YZhang Y HShao M2006Atmospheric Environmental ChemistryBeijingHigher Education Press566756–67(in Chinese)
2Yao LLiu W QLiu J GKan R FXu Z YRuan JDai Y H 2015 Chin. J. Lasers 42 0215003 (in Chinese)
3Logan J APrather M JWofsy S CMcelroy 1981 Journal of Geo-physical Research 86 7210
4Change C2007Synthesis Report. Contribution of Working Groups I, II and III to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change [J]. IPCC 2007
5http://www.wmo.int/pages/prog/arep/gaw/documents/Final_GWA_213_web.pdf
6Hu Y Y2012Physics41495
7Willer USaraji MKhorsandi AGeiser PSchade W 2006 Optics and Lasers in Engineering 44 699
8Pan W DZhang J WDai J MSong K2012Spectroscopy and Spectral Analysis3210
9Kan R FLiu W QZhang Y JLiu J GDong F ZWang MGao S HChen JWang X M2005Chin. J. Lasers321217(in Chinese)
10Chao XJeffries J BHanson R K 2009 Measurement Science and Technology 20 115201
11Xia HDong F ZWu BZhang Z RPang TSun P SCui XHan LWang Y 2015 Chin. Phys. 24 034204
12Gabrysch MCorsi CPavone F SInguscio M 1997 Appl. Phys. 65 75
13Engelbrecht R 2004 Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy 60 3291
14Chen DLiu W QZhang Y JLiu J GWei Q NKan R FWang MCui Y BChen J 2006 Chin. Phys. Lett. 23 2446
15Wang T DLiu W QZhang Y JWang M2007Infrared and Laser Engineering36336(in Chinese)
16Tang Y YLiu W QKan R FZhang Y JChen DZhang SRuan J 2012 Chin. Opt. Lett. 10 041404
17Yan F LZhang J CYao D YLiu F QWang L JLiu J QWang Z G 2015 Chin. Phys. 24 024212
18Moeskops B WCristescu S MHarren F J 2006 Opt. Lett. 31 823
19Kosterev AWysocki GBakhirkin Yet al. 2008 Appl. Phys. 90 165
20Wysocki GBakhirkin YSo STittle F KHill R QFraser M P 2007 Appl. Opt. 46 8202
21Mappé IJoly LDurry GThomas XDecarpenterie TCousin JDumelie NRoth EChair AGrillon P G 2013 Review of Scien-tific Instruments 84 023103
22Tao LSun KKhan M AMiller D JZondlo M A 2012 Opt. Express 20 28106
23Cao Y CSanchez N PJiang W ZGriffin R JXie FHughes L CZah CTillel F K 2015 Opt. Express 23 2121
24Li HRieker G BLiu XJeffries J BHanson R K 2006 Appl. Opt. 45 1052
25Klein AWitzel OEbert V 2014 Rapid Sensors 14 21497
26Fernholz TTeichert HEbert V 2002 Appl. Phys. 75 229
27Zhang Z RWu BXia HPang TWang G XSun P SDong F ZWang Y 2013 Acta Phys. Sin. 62 234204 (in Chinese)
28Rothman L SGordon I EBabikov Yet al. 2013 Journal of Quantitative Spectroscopy and Radiative Transfer 130 4